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1BSTBACT 


Horizontal  and  vortical  theraal  structures  were  examined 
in  a  region  of  the  northeastern  Pacific  Ocean  off  Northern 
California.  The  observations  were  acquired  on  three  cruises 
as  part  ,of  the  ONfi-sponsored  Ocean  Prediction  Through 
Observation,  Modeling  and  Analysis  (OPTOMA)  research 
progcaa#  (centered  in  a  region  of  the  California  Current 
System,  ca.  37  to  39«B,  124  to  126°N,  during  June  and  July, 
1983 . 

The  horizontal  temperature  correlation  scale  was  between 
30  and  50km,  which  was  a  significant  factor  when  comparisons 
were  made  between  measured  horizontal  theraal  structures  and 
those  retrieved  from  Fleet  Numerical  Oceanography  Center's 
(FNOC)  analyses,  which  had  grid  spacing  of  320km  (TEOTS)  and 
40ka  (EOIS) . 

Operational  (modeled)  analysis  fields  were  not  in  agree¬ 
ment  with  the  obsexved  fields.  The  major  discrepancies 
occurred  in  the  magnitude  of  the  sized  layer  depth  and  the 
shape  of  the  horizontal  temperature  fields  (maps) . 

The  operational  significance  of  differences  between 
modeled  and  observed  thermal  structures  was  assessed  in 
terms  of  their  effect  on  low-frequency;  i.e.,less  than  IKhz, 
acoustic  propagation  utilizing  the  FACT9H  and  PE 
Transmission  Loss  models. 

Mixed  layer  depth  differences  produced  significant  disa¬ 
greements  between  direct  ranges  computed  from  model  and 
observed  temperature  profiles.  The  effect  was  most 
pronounced  at  higher  frequencies  and  when  both  source  and 
receiver  were  shallow;  i.e. ,  both  at  20m. 

A  comparison  was  made  between  average  depth/temperature 
profiles  from  July,  1982,  July,  1983,  and  FNOC  climatology 


to  obtain  a  ae&sure  of  the  effect  of  interannual  variability 
in  the  doaain.  This  coaparison  shoved  that  a  significant 
tespexature  anosaly  existed  in  the  upper  400s  in  1983 
coapaxed  to  1982  due  to  El  Sino  and  that  this  anoaaly  vas 
not  xepresented  by  tfce  FNOC  cliaatology. 

The  differences  between  aodeled  and  aeasured  theraal 
structures  axe  believed  to  be  related  to  theraal  structure 
aodel  resolution,  aodel  sensitivity  to  input  data,  short 
scales  of  spatial  variability  and  non- representative  cliaa¬ 
tology  for  the  doaain. 
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I*  I1TBODOCTIOI 

1.  B1CKGHOOID  AID.  OBJECTIVES 

The  dependence  of  acoustic  propagation  on  the  ocean 
theraal  structure  and  the  increased  deaand  fro*  the  Navy  for 
accurate  acoustic  forecasts  to  facilitate  planning  and 
econoaical  use  of  resources  has  led  to  an  increasing 
reguireaent  for  a  reliable  and  accurate  aodel  to  represent 
ocean  theraal  structure.  Such  a  reguireaent  is  not  easily 
net,  however,  because  of  the  United  nuaber  of  observations 
in  the  typical  ocean  theraal  structure  data  base,  approxi- 
aately  150  XBT  measurements  per  day  in  the  northern  heai- 
sphere  [Clancy  et  al.,1983].  In  view  of  the  upper  ocean's 
spatial  and  tenporal  scales  of  variability  such  a  snail  data 
base  represents  a  significant  handicap  to  the  perforaance  of 
any  ocean  analysis  acdel. 

The  availability  of  an  ocean  theraal  structure  data  base 
froa  a  relatively  densely  saapled  area  of  the  Northeastern 
Pacific  Ocean  (approxinately  200ka  vest  of  point  Arena, 
Northern  California)  aade  it  possible  to  evaluate  the  accu¬ 
racy  of  the  output  of  two  operational  ocean  theraal  struc¬ 
ture  analysis  systens. 

-The  objectives  of  this  research  are  firstly  to  obtain  a 
description  of  the  characteristics  of  the  ocean  therxal 
structure  froa  a  portion  of  the  ocean  (ca.37°N  to  39«N; 
12h°H  to  1260W)  saapled  during  the  three  cruises  conprising 
CPT0MA5  which  took  place  in  June  and  July  1983.  (OPTOMA  - 
Ocean  Prediction  Through  Observations,  Modeling,  and 
Analysis,  is  a  joint  Harvard/NPS  project  'intended  to 
acguire  field  data  to  characterise  synoptic  scale  eddies 
over  a  donain  in  the  California  Current  off  Northern 
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California,  and  to  'set-up*  an  eddy-resolving,  statistical/ 
dynamical,  United  domain,  open  boundary  nuaerical  ocean 
prediction  model. 'Tie  *5*  indicates  the  fifth  series  of 
cruises  of  the  project).  Secondly,  the  analysed  ocean 
thermal  structure  is  to  be  used  as  a  'sea  truth'  by  which  to 
evaluate  the  accuracy  of  two  real-time  ocean  thermal  anal¬ 
ysis  systems  presently  in  use  at  the  Fleet  Numerical 
Oceanography  Center  (FNOC) ,  namely,  the  Expanded  Ocean 
Thermal  Structure  (FCTS)  analysis  and  the  Thermodynaiic 
Ocean  Prediction  System  (TOPS) -Coupled  EOTS  (TEOTS)  anal¬ 
ysis.  The  final  objective  is  to  determine  the  operational 
significance  of  deviations  in  the  modeled  thermal  structure 
from  the  'sea  truth'  by  utilizing  the  output  from  opera¬ 
tional  acoustic  models  to  provide  a  quantitative  measure  of 
the  effect  of  these  thermal  structure  differences  on 
acoustic  propagation  in  the  ocean.  Two  low-f reguency, 
passive-acoustic  models  were  used  to  provide  a  measure  of 
the  variability  of  acoustic  energy  propagation  which  could 
be  expected  from  the  acoustic  models  alone.'  These  then, 
would  provide  a  benchmark  by  which  the  variability  due  to 
differences  between  modeled  and  analysed  thermal  structures 
could  ie  judged  significant  or  not.  The  two  acoustic  models 
used  were  the  FACT9H  (Fast  Asymptotic  Ccherrent 
Transmission)  transmission  loss  model  and  the  AESD  (Acoustic 
Enviromental  Support  Detachment)  PE  (Parabolic  Eguaticn) 
model  [  Brock,  1 978  ]. 

B.  STUDY  ABBA  AID  ITS  OCEANOGBAPHI 


The  area  investigated  is  located  in  the  California 
Current  System  (CCS) ,  south  of  the  Mendocino  Escarpment.  The 
region  is  approximately  bounded  by  37®  to  39°  N  and  124°  to 
126°N.  The  bathymetry.  Fig.  1.1  ,  shows  the  area  to  be 
offshore  of  the  continental  slope,  with  an  average  water 
depth  of  4000m. 


Figure  1.1  Bathyaetry  of  Study  Doaain  and  Adjacent  Area 


the  California 


The  CCS  is  coaposed  of  four  currents: 

Current,  the  Davidson  Current,  the  California  Undercurrent, 
and  the  Southern  California  Current.  Of  these  four  currents 
the  southward  flowing  California  Current  is  the  aajor  influ¬ 
ence  in  the  study  dcnain.  This  current  forns  the  eastern 
boundary  current  of  the  large  anticyclonic  gyre  centred  near 
the  Hawiian  Islands,  As  an  extension  of  the  Nestwind  Drift, 
the  California  Current  occurs  between  the  North  Pacific 
ataospheric  high  pressure  systea  and  the  seai-peraanent 
theraal  low  positioned  over  Central  California.  Contrary  to 
the  classical  picture  of  a  broad,  slow,  shallow  cliaatolo- 
gical  aean  eastern  boundary  current,  the  CCS  is  known  to 
coaprise  aesoscale  eddies,  aeandering  jets  and  turbulent 
filaaents  [Hooers  and  Robinson,  1984]  together  with  their 
associated  coaplex  horizontal  and  vertical  theraal  struc¬ 
tures.  The  aean  current  is  broad  (approx.  lOOOka  wide), 
shallow  (less  than  500a  deep)  and  sluggish  aoving  (approx 
0.9ka/hr).  The  southward  to  southeastward  winds,  which 
prevail  during  the  spring  to  the  fall  months  south  of  40®  S, 
produce  upwelling  along  the  coast,  extending  50  to  100ka 
offshore.  Upwelling  begins  in  February  or  Rarch  and 
continues  through  August  in  the  study  doaain.  The  cold, 
usually  high  salinity  water,  brought  into  the  surface  layers 
by  this  process,  is  then  aoved  offshore  by  the  saaller  scale 
perterbations  in  the  CCS,  further  coaplicating  the  ocean 
theraal  structure.  The  aost  intense  upwelling  sites  occur  to 
the  south  of  capes  and  headlands  [Reid  et  al.,1958]  with  the 
aost  conspicuous  centres  of  upwelling  adjacent  to  the  study 
doaain  at  35°  N  and  41®  N  [Sverdrup  et  al.,  1942  and  Reid  et 
al. ,  1958]. 

Mean  aonthly  upwelling  velocities  for  the  west  coast  of 
North  AaericafFig  1.2;  froa  Heath  (1983)),  based  on  Nelson's 
(1977)  wind  stress  curl  data,  clearly  show  upwelling  to  be 
at  its  aost  intense  in  July  off  Northern  California.  Hence, 


the  aaxiaua  horizontal  and,  because  of  the  differing  water 
■asses  involved,  vertical  variations  in  theraal  structure 
should  be  expected  to  occur  in  the  study  doaain  during  the 
period  under  investigation. 

In  addition  to  the  seasonal  variations  in  theraal  struc¬ 
ture,  inter-annual  variability  is  also  evident.  Significant 
increases  in  temperature  off  Central  and  Northern  California 
occurred  during  the  winters  of  1972,  1976,  1979  and  1983, 
corresponding  to  El  Nino  episodes  [Breaker,  1983].  Of 
particular  interest  are  the  effects  of  the  aajor  El  Nino 
waraing  of  1982-83.  1  brief  analysis  of  the  possible  opera¬ 
tional  ccnseguences  of  such  an  anoaolous  event  is  presented 
in  Chapter  V. 


UPWELLING  VELOCITY  {IO'3cm/sec) 


Figure  1.2  lean  Honthly  Opwelling  Indices. 
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II.  DBSCBIPTIOH  OJ  COMPOTES  MODELS 
A.  OC1A1  ABALYSIS  MODELS 

1.  Expanded  Ocean  Thermal  Structure  (SOTS)  Analysis 

Model 

The  Expanded  Ccean  Thermal  Structure  (EOTS)  analysis 
[Mendenhall  et  al.,  1978]  has  been  the  OS  Navy’s  operational 
thermal  analysis  system  for  the  past  several  years  and 
provides  ocean  thermal  structure  input  to  most  of  the 
acoustic  predictions  generated  at  FNOC. 

EOTS  is  run  cn  a  63x63  grid  for  both  hemispheric  and 
regional  coverage.  The  grid  spacing  for  the  hemispheric 
product  is  200  to  400km,  while  that  of  the  regional  version 
is  typically  40  to  80kn.  For  the  SOCAL  region,  in  which  the 
study  domain  resides,  the  spacings  are  320km  and  40km 
respectively. 

The  EOTS  analysis  system  is  based  on  a 
Tields-by- Information-Blending  (FIB)  methodology  [Boll  and 
Mendenhall,  1971;  Hell,  Cuming  and  Mendenhall,  1979]  and  is 
a  numerical  blending  technigue  whereby  the  difference 
between  a  first  guess  field  (in  this  case  the  previous  anal¬ 
ysis  weighted  towards  climatology  at  a  specied  rate)  and  the 
assembly  field,  are  eliminated  in  a  weighted,  least  mean 
squares  sense.  The  first  guess  field  provides  shape  informa¬ 
tion  for  the  horizontal  blending  procedure.  The  assembly 
field,  into  which  all  new  observations  are  assembled,  is 
obtained  in  a  two-step  process.  First,  the  assembly  field 
from  the  previous  analysis  is  trended  towards  climatology  in 
the  same  manner  as  tte  first  guess  field,  and  a  weight  field 
associated  with  it  is  decayed  an  arbitrary  amount.  Then,  a 
weighted  average  of  the  assembly  and  first  guess  fields  is 


FIXED  LEVELS 


DEPTH  (m) 


The  vertical  grid  on  which  the  twenty-six  EOTS 
parameters  are  defined.  Columns  labeled  T,  T',  T"  are 
temperature,  first  vertical  temperature  derivative,  and 
second  vertical  temperature  derivative,  respectively. 
Parameters  y-2b  are  associated  with  fixed  levels. 
Parameters  2-8  are  associated  with  floating  levels 
defined  relative  to  parameter  1,  Primary  Layer  Depth 
(PLD).  _ 


Figure  2.1  Analysis  Levels  for  FNOC  Analyses. 
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formed  (using  a  low,  specified  weight  for  the  first  guess 
field)  to  provide  a  new  assembly  field.  New  observations  are 
interpolated  to  the  nearest  analysis  grid  point  using  shape 
inforsation  froa  the  first  guess  field.  At  each  grid  point, 
a  weighted  average  of  the  interpolated  observations  and  the 
asseably  field  is  formed  and  a  final  assembly  field  deter¬ 
mined.  It  is  upon  this  final  asseably  field  and  the  first 
guess  field  that  the  FIB  analysis  is  performed.  The  proce¬ 
dure  is  designed  to  spread  information  from  data-rich  points 
to  data-sparse  points  and  produces  a  smoothed,  horizontally 
and  vertically  blended  field.  During  the  blending,  both 
first  guess  and  assembly  fields  are  weighted  towards  clima¬ 
tology,  and,  in  the  absence  of  inputs,  they  will  revert  to 
climatology  over  a  period  of  approximately  fifteen  days 
[Pollack,  personal  communication]. 

The  analysis  is  performed  on  twenty-six  parameters 
at  fixed  and  floating  levels  from  the  surface  to  400a,  Fig 
2. 1  .  Parameter  1  is  the  Primary  Layer  Depth  (PLD)  which  in 
general  equates  to  the  depth  of  the  surface  mixed  layer.  The 
remaining  parameters  are  temperature  at  selected  depths  and 
vertical  temperature  gradients.  Parameters  2  to  8  are  anal¬ 
ysed  at  floating  levels  defined  relative  to  the  PLD. 
Parameters  9  to  26  are  associated  with  fixed  levels.  Below 
400m  the  thermal  field  is  constructed  from  climatological 
archives  modified  to  blend  smoothly  with  the  analysed  temp¬ 
erature  profiles  above. 

Selection  of  the  PLD  (point  1)  is  accomplished  by 
applying  an  algorithm  to  all  BT  observations  made  during  a 
specified  number  of  days  (typically  five)  prior  to  the  anal¬ 
ysis  time  to  identify  PLD  candidates  (i.e. ,  the  depth  (s)  at 
which  the  criteria  prescribed  by  the  PLD  algorithm  are 
satisfied)  from  each  profile.  The  algorithm  assigns  the  PLD 
to  the  depth  (s)  (below  the  surface)  of  maximum  rate  of 
curvature  in  the  vertical  profile.  Then  a  preliminary 


three-cycle  FIB  is  performed  using  the  selected  FID 
candidates  as  input,  and  the  PLD  and  assembly  fields  from 
the  previous  analysis,  trended  towards  climatology,  as  the 
first  guess  and  initial  assembly  fields.  Any  recent  BT's  are 
then  examined  and  the  PLD  candidate  from  each  that  is 
closest  to  the  PLD  produced  by  the  preliminary  FIB  analysis 
is  selected  and  entered  into  a  second  three-cycle  FIB 
analysis.  This  produces  the  final  analysed  PLD. 

2.  TOPS-Coupled  EOTS  (TEOTS)  Analysis 

The  TEOTS  system  is  an  attempt  to  incorporate 
oceanic  physics  and  air-sea  interactive  processes  into  the 
analysis  by  coupling  the  ocean  thermal  analysis  to  the 
atmospheric  forcing  via  the  physics  of  the  TOPS  mixed  layer 
model. 

The  TOPS  mixed  layer  model  is  the  Navy's  vertical 
ocean  thermal  structure  forecast  system.  Given  that  the 
ocean  is  primarily  atmospherically  forced,  and  that  mixed 
layer  depths  have  proved  'highly  predictable  with  a  variety 
of  models',  such  as  the  Denman  (1973)  or  Mellor  and  Durbin 
(1975)  models,  TOPS  was  developed  by  NORDA  (Naval  Ocean 
Research  and  Development  Activity) ;  it  is  based  largely  on 
the  use  of  rather  well-defined  atmospheric  variables  to 
drive  a  model  to  describe  the  relatively  data  sparse  oceans. 
Surface  wind,  solar  heat  fluxes,  and  precipitation  fields 
are  among  those  necessary  to  drive  the  TOPS/TEOTS  models, 
i  The  surface  fluxes  are  supplied  by  NOGAPS  (Naval  Operational 

Global  Atmospheric  Prediction  System)  and  the  physics  of  the 
model  are  based  on  the  Mellor  and  Yamada  (1974)  level-2 
turbulence  model  with  the  added  effects  of  advection  by 
instantaneous  wind  drift  and  climatologically  averaged  geos- 
trophic  currents.  The  TOPS  model  produces  a  twenty-four 
hour  forecast  for  temperature  fields  at  seventeen  depths 
from  the  surface  to  500m  on  a  horizontal  grid  spacing 
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identical  to  that  used  in  the  teots  analysis  to  which  i 
coupled.  The  main  purpose  behind  the  TOPS/TEOTS  approac 
to  de-couple  the  objective  analysis  component  (EOTS)  SSI 
MID  fields  from  climatology,  and  to  provide  realistic  ] 
boundary  conditions  fcr  the  atmospheric  models. 

TEOTS  is  coupled  to  TOPS  in  cyclical  fas] 
providing  initial  conditions  for  a  24hr  TOPS  forecast 
is  subsequently  fed  back  into  EOTS  as  a  first  guess  i 
for  the  following  day's  analysis.  With  the  exception  oi 
lock-step  procedure  described  above  and  a  diff« 
prescription  for  certain  parameters,  the  TEOTS  ana] 
system  is  the  same  as  EOTS,  relying  on  the  FIB  technigi 
combine  the  various  types  of  data.  Clancy  and  Pollack  (‘ 
give  a  mere  rigorous  description  of  both  EOTS  and  TEOTS, 

For  the  study  domain  EOTS  was  used  for  the  reg; 
(NOCAL)  analysis  and  TEOTS  for  the  hemispherical  analys] 

B.  ACOUSTIC  MODELS 

1.  The  FACT9H  Model 

The  FACT9H  Transmission  Loss  Model  is  the 
Interim  Standard  Transmission  Model  for  ocean  regions  < 
may  be  treated  with  a  single  sound-speed  profile  and  a 
ocean  bottom  [Jacobs,  1982].  It  is  a  ray-acoustics  i 
designed  for  the  computation  of  transmission  loss  as  a 
tion  cf  range  and  frequency  for  a  source  and  receiv 
fixed  depths.  The  classical  ray  treatment  has  been  augm 
with  the  addition  of  higher  order  asymptotic  correctio 
the  vicinity  of  caustics,  and  the  phase  addition  of  ce 
ray  paths. 

Wave  theory  for  acoustic  propagation  states 
there  is  a  frequency  below  which  no  'trapped'  normal 
may  exist  for  a  given  channel  or  duct  (waveguide)  [Cc 
et  al.,  1980].  The  FACT9H  model  attempts  to  translate 
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condition  to  ray-acoustics  by  applying  an  algorithm  [ Clay, 
1968]  to  characterize  the  gross  features  of  ducted  propaga¬ 
tion.  The  intensity  in  the  surface  duct  is  found  from  the 
principle  of  conservation  of  energy  modified  by  additional 
losses  (proportional  to  range)  caused  by  duct  leakage  and 
rough-surface  scattering  of  energy  from  the  duct  [Marsh  and 
Schulkin,  1967].  For  both  source  and  receiver  within  the 
surface  duct  the  transmission  loss  (TL)  as  a  function  of 
range  (R)  is  given  by  Egn.  2.1 

TL  (R)  =22  ♦  lOlog  (R  .0)  ♦  b.R  (egn  2.1) 

where  8  is  the  angle  for  the  ray  at  the  surface  which  just 
grazes  the  bottom  of  the  duct  and  b  contains  the  duct- 
leakage  and  rough-surface  losses, Egn.  2.2 

b=  1 4. 88  x  105  f-s/3  g-i/3  z-3  ♦  *  £1/2  z-yz  (egn  2.2) 

In  the  above  expression  f  is  the  freguency,  g  the  magnitude 
of  the  below  layer  gradient,  z  the  duct  depth,  and  w  is  a 
factor  determined  by  surface  wave  height. 

The  9H  is  a  model  designator  and  reflects  the  stage 
of  development  of  the  model  physics  (FACT  10  is  presently  in 
operation  at  FNOC) . 

2 .  The  PE  Model 

The  PE  Model  replaces  the  reduced  elliptic  wave 
equation  with  a  paralclic  partial  differential  eguation  that 
can  be  integrated  numerically  using  the  Tappert-Hardin 
split-step  Fourier  Algorithm  [Brock,  1978].  The  parabolic 
wave  eguation  includes  diffraction  and  all  other  full-wave 
effects  as  well  as  range  depedent  environments.  The  entire 
range  and  depth  dependent  acoustic  field  is  computed  as  the 
solution  is  marched  forward  in  range.  The  model  assumes  a 
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flat  pressure  release  ocean  surface  and  a  vanishing  field  at 
the  depth  of  the  finite  Fourier  transform  (i.e.,  for  bottom 
grazing  angles  greater  than  the  maximum  permissible  33° ,  the 
acoustic  energy  is  smoothly  attenuated  to  zero  and  back 
reflection  of  acoustic  energy  into  the  water  column  is 
avoided) .  For  grazing  angles  less  than  the  maximum,  a 
pseudo  radiation  condition  is  introduced  at  the  water-bottom 
interface  by  smoothly  attenuating  the  field.  Since  the  error 
in  the  parabolic  approximation  increases  as  angles  increase 
from  the  horizontal,  steep  bottom  slopes  can  cause  inaccura¬ 
cies.  The  model  is  primarily  considered  useful  for 
predicting  low-freguency  acoustic  propagation  of  energy 
along  waterborne  or  shallow-angle,  bottom-bounce  paths  in 
deep  water. 


III.  fill  £152 


1.  DATA  ACQUISITION 

The  XBT  and  CTD  data  from  the  surface  to  500a  acquired 
froa  the  OPTOMA5  cruises  foraed  the  primary  source  of  field 
data.  CPTCMA5  consisted  of  three  separate  cruises  undertaken 
between  15  June  and  19  July,  1983  by  R/  V  ACANIA.  The  three 
cruises,  designated  AI,  All  and  AIII  each  saapled  approxi- 
aately  the  same  area  of  the  CCS,  Fig.  3.1 

A  secondary  source  of  field  data  was  provided  by  deep 
CTD  casts  froa  the  cruise  undertaken  by  the  R/7  DE  STEIGUER 
between  10  and  21  June,  1983.  A  section  of  this  cruise 
saapled  the  OPT OH A5  area  and  provided  CTD  data  to  3000a. 
Satellite  imagery  was  also  used,  when  available,  to  provide 
a  qualitative  measure  of  the  sea  surface  temperature 
distribution. 

Temperature  profiles,  to  be  compared  later  with  the 
field  data,  were  retrieved  from  the  archives  of  FNOC  for  the 
same  area  and  period.  Both  EOTS  and  TEOTS  depth /temperature 
profiles  for  the  beginning  and  end  of  each  OPTOMA5  cruise 
were  obtained  for  selected  positions.  Fig.  3.2  . 

Details  of  the  data  base  are  given  in  Table  I  . 

B.  DESCRIPTION  OF  DATA  AID  DATA  PROCESSING 

1.  ACANIA  £a£a 

The  positions  of  the  XBT/CTD  casts  for  the  Acania 
cruises  are  documented.  Figs.  3.3  to  3.5.  All  XBT/CTD  data 
were  edited  thus  eliminating  any  obviously  erroneous 
profiles.  A  secondary  visual  inspection  of  each  profile  was 
undertaken  with  the  result  that  a  final  data  base  for  each 
cruise  was  obtained,  Tablo  i 


TABLE  I 

Details 

of  Data  Base 

SOURCE 

PERIOD 

MO.  OF  ACCEPTABLE 
PROFILES 

RANGE  OF 
PROFILES 

OPTOMA  AI 

6/15-6/22 

74 

0-5Q0m 

OPTOMA  All 

6/29-7/4 

89 

0-50  0m 

OPTOMA  AIII 

7/13-7/19 

86 

0-50  0m 

DI  STEIGUER 

6/10-6/21 

8 

0-30  00  m 

EOTS/TEOTS 

6/16,6/23, 
6/30*7/7. 
7/ 10*7/21 

42 

0-40 0m 

The  XBT's 

were  then 

interpolated  in  the 

vertical  to 

standard  5b  depths  fron  the  surface  to  500m.  The  loss  of 
teaperature  microstructure  resulting  from  this  interpolation 
uas  assumed  unimportant  in  view  of  the  goals  of  this  study. 
These  edited,  *  smoothed*  depth/temperature  profiles  formed 
the  field  data  upon  which  the  later  analysis  was  to  be 
based. 

Information  from  each  cruise  was  assumed  'quasi- 
synoptic*  (i. e.  the  region  was  believed  to  have  varied  some¬ 
what  over  the  period  of  sampling) .  This  assumption  is 
weakest  in  the  near  surface  layers  where  the  time  scale  of 
variability  can  be  of  the  order  of  a  day  or  less.  The 
validity  of  this  assumption  in  the  context  of  the  present 
study  is  addressed  in  Chapter  V. 

2.  fi  STEIGUEB  fia^a 

The  DE  STEIGDER  sampling  provided  CTD  casts  from  the 
surface  to  a  maximum  depth  of  3000m.  These  profiles  were 
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also  interpolated  vertically  to  5»  depth  intervals  and, 
based  on  the  assumption  that  the  water  column  was  near 
stationary  below  the  thermocline  (approximately  500m) ,  they 
were  used  to  extend  the  ACANIA  and  FNOC  profiles  to  a  depth 
of  3000m.  These  profiles  were  then  extrapolated  linearly  to 
the  ocean  bottom  (assumed  flat  at  4000m) . 

An  average  salinity/depth  profile  was  calculated 
from  the  ’in-area'  El  STEIGUER  CTD  casts  and  this  average 
profile  was  assumed  as  the  standard  salinity  profile  for  all 
sound  speed  calculations. 

3 .  Satellite  I magerv 

Advanced  Very  High-Resolution  Radiometer  ( AVHRR) 
infrared  images  from  the  NOAA-7  satellite  were  used  to  qual¬ 
itatively  show  the  mere  obvious  features  of  the  SST  field. 
These  images  were  obtained  on  an  'as  available'  basis  during 
each  of  the  ACANIA  cruises.  Details  of  the  major  features 
shown  in  these  images  were  sketched.  Figs.  3.6  to  3.8  The 
satellite  images  were  also  used  to  gualitively  assess  the 
assumption  of  quasi-synopticity. 

4.  FNOC  Archived  Data 

Temperature  profiles  for  the  standard  positions. 
Fig.  3.2  ,  were  extracted  from  the  archives  of  FNOC.  These 
profiles  were  used  for  both  EOTS  and  TEOTS  analyses  on  the 
dates  closest  to  the  commencement  and  termination  of  each 
ACANIA  cruise. 


Pigare  3.6  Satellite  Observed  SST  Features-Craise  il 


Figure  3.7  Satellite  Observed 


Figure  3.8  Satellite  Observed  SST  Features-Cruise  AII1 
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To  facilitate  comparison  with  FNOC's  data  fields, 
study  domain  was  divided  into  a  7x7,  sguare  grid  w: 
grid-spacing  of  20km  and  the  irregularly  spaced  XB' 
profiles  interpolated  onto  this  grid. 

The  horizontal  correlation  of  the  temperature  fie; 
various  depth  intervals  in  the  top  150m  of  the  water  c< 
were  evaluated.  The  correlation  at  various  levels  disp; 
similar  chracteristics  showing  zero  crossings  occurrii 
approximately  80km  with  a  0.5  correlation  at  approxim< 
30km  (approximately  the  e-folding  distance  between 
origin  and  the  zero  crossing  distance).  It  was,  there; 
reasonable  to  assume  that  the  temperature  variations 
points  within  30  km  cf  a  given  position  within  the  d< 
were  related  (i. e.  all  points  within  30km  could  be  exp< 
to  influence  a  given  point) .  As  an  example,  the  SST  ] 
zontal  correlations  for  each  cruise  are  presented.  Figs, 
to  4.3  .  The  correlation  information  thus  derived  was 
when  interpolating  the  field  data  onto  the  grid, 
observations  within  30km  of  each  grid  point  were  usei 
each  observation  value  was  weighted  according  to  the  in1 
of  its  distance  from  the  grid  point.  The  interpolated 
thus  formed  were  then  automatically  contoured  at  FNOC  i 
erature  levels.  Fig.  2.1  The  mixed  layer  depth  (HLD)  a 
lated  from  the  field  data  was  equated  to  the  PLD  in  the 
analyses  (point  1  Fig.  2.1)  and  was  defined  as  the  s 
lowest  depth  below  which  a  negative  temperature  gradieni 
found  with  a  magnitude  greater  than  3°C  /100m  over 

consecutive  5m  depth  intervals.  The  algorithm  usee 


calculate  the  MLD  was  designed  to  yeild  the  true  MLD  even  in 
the  presence  of  strong,  near-surface  negative  temperature 
gradients  caused  by  diurnal  heating  the  "afternoon  effect" 
which  seldom  exceeded  7m  in  depth.  MLD's  calculated  by  this 
method  agreed  well  with  manually  determined  MLD’s  for  the 
sampled  profiles.  FNCC  EOTS  and  TE01S  profiles  were  also 
interpolated  and  contoured  in  a  similar  manner  and  on  the 
same  grid.  The  HIE  computed  for  the  model  profiles  using 
the  algorithm  described  above  was  in  agreement  with  the  PID. 

By  contouring  the  fields  in  the  manner  described,  it  was 
possible  to  extract  a  depth/temperature  profile  from  the 
surface  to  4000m  for  any  position  within  the  study  domain. 

Throughout  the  study  domain,  temperature  variations  at 
400m  were  small  enough  (<0.5°C)  for  all  profiles  (including 
those  provided  by  FMOC)  to  be  smoothly  blended  between 
shallow  and  deep  profiles  above  the  depth  of  the  deep  sound 
channel  axis  (ca.  500m)  . 

The  feasibility  of  using  discrete  depth/temperature 
levels,  with  depth  separations  dictated  by  the  EOTS  analysis 
levels,  to  represent  the  vertical  temperature  structure  of 
the  area  was  tested  by  examining  the  vertical  pattern  corre¬ 
lation  between  temperatures  in  the  near  surface  layers  (the 
zone  of  maximum  variation  of  temperature  with  depth) ,  Figs. 
4.4  to  4.6  .  These  pattern  correlations  are  similar  for  each 
cruise  and  have  zerc  crossings  at  depths  of  the  order  of 
300m  when  the  SST  is  correlated  with  deeper  levels.  Such  a 
pattern  correlation  gives  an  good  indication  of  the  prob¬ 
ability  that  variations  in  temperature  at  one  depth  are 
related  to  variations  in  temperature  at  another  depth.  At 
the  e-folding  depth  (approximately  100m) ,  the  correlation  is 
approximately  0.7.  It  was  therefore  considered  reasonable  to 
represent  the  vertical  temperature  structure  at  the  discrete 
levels  dictated  by  the  EOTS  analysis.  The  EOTS  depth 
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Figure  4.2  SST  Horizontal  Correlation~Cruise  ill. 
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layer  is  of  the  order  of  50b  which  yields  a  pattern  correla¬ 
tion  coefficient  ranging  froa  0.6  for  cruise  All  to  0.8  for 
cruises  AI  and  AIII.  Belov  the  nixed  layer  (approxiaately 
50a  )  the  correlation  improves  to  0.8  or  better  between 
adjacent  contoured  levels. 

The  smoothing  effect  of  deriving  profiles  as  described 
above  (i.e.  the  higher  frequency  fluctuations  of  tenperature 
with  depth  are  reaoved)  was  not  considered  important, 
because  the  profiles  thus  derived  were  later  to  be  used  as 
inputs  tc  low  frequency  (<1000Hz)  acoustic  nodels  and  the 
essential  shape  cf  the  profile  was  maintained. 
Additionally,  ray  theory  is  only  valid  if  sound  velocity 
does  not  change  much  in  a  wavelength,  therefore,  the  inclu¬ 
sion  of  the  higher  wavenuaber  fluctuations  of  temperature 
(and,  therefore,  sound  speed)  with  depth  could  invalidate 
the  use  of  this  theory  [Brick,  1983  ]. 

Four  features  of  the  vertical  temperature  profile  were 
selected  as  being  of  primary  importance  in  describing  the 
thermal  structure  for  low  frequency  passive  acoustic  appli¬ 
cations.  These  were:  sea  surface  temperature  (SST) ,  mixed 
layer  depth  (MLD)  (previously  defined) ,  the  below  layer 
gradient  (BLG)  (defined  as  the  temperature  gradient  in  the 
25m  immediately  below  the  MID)  and  the  thermocline  gradient 
(THG)  (Defined  as  the  temperature  gradient  between  the  MLD 
and  the  deep  sound  channel  axis  depth) .  The  contour  fields 
derived  for  each  of  these  parameters.  Figs.  4.7  to  4.10  , 
are  described  in  detail  in  Section  B  of  this  chapter. 

B.  DESCBIPTIOH  OF  ABE1  THERMAL  STRUCTURE 

The  following  is  a  brief  description  of  the  thermal 
structures  measured  in  the  study  domain  during  the  three 
CPT0MA5  cruises.  The  desciption  in  limited  to  the  four 
parameters  defined  in  section  A  above. 
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Fiance  4.7  SST  Contours 
a.  Cruise  AI  b.  Cruise  III  c.  Cruise  A  III. 
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The  significant  weather  for  each  cruise  (i.e.  wind  speed 
and  cloud  cover)  were  obtained  from  R/V  ACANIA’s  hourly  log. 
The  wind  was  from  a  generally  southeastward  throughout  all 
three  cruises,  and,  for  the  aost  part,  wind  speeds  were 
between  15  and  20  knots.  The  major  excursions  from  this  norm 
were  (a)  at  the  beginning  of  cruise  AI  when  wind  speeds  were 
greater  than  30  knots  on  the  first  and  second  days  (16  and 
17  Jun),  (b)  during  the  latter  part  of  cruise  AI  and  early 
in  cruise  All  when  winds  of  10  to  15  knots  were  recorded  and 
(c)  during  the  first  two  days  of  cruise  AIII  when  wind  speeds 
again  increased  to  30  knots  or  more.  Throughout  the  three 
cruises,  partly  clcudy  to  clear  skies  prevailed,  with 
predominantly  less  than  4  oktas  of  cloud  cover.  Pressure 
remained  nearly  steady  at  approximately  1030  mb. 

1.  Cruise  A I  (1.5  Jun-22  Jun,  1  983) 

SST’s  ranged  from  14.60C  in  the  northwest  of  the 
domain  tc  12.2°C  associated  with  a  cool  feature  protruding 
into  the  southeast.  Fig.  4.7  (a).  The  effect  of  the  cool 
intrusion  was  to  intensify  the  horizontal  temperature 
gradient  in  its  vicinity  and  to,  locally,  distort  the  other¬ 
wise  nearly  meridional  SST  structure.  The  maximum  tempera¬ 
ture  gradient  intensification  occurred  on  the  northern  flack 
of  the  cool  feature,  where  a  gradient  of  7.0°C/100m  was 
observed.  The  position  of  the  cooler  water  agreed  well  with 
the  position  obtained  from  satellite  imagery.  Fig  3.6. 

The  depth  of  the  mixed  layer.  Fig.  4.8  (a)  ,  varied 
from  20  to  55m  in  a  complex  manner  throughout  the  area  with 
the  shallowest  depths  in  a  50km-wide  swath  southeast  to 
northwest  across  the  centre  of  the  area.  No  discernatle 
correlation  was  observed  between  the  SST  and  MLD  patterns. 

Below  layer  gradients  ranged  from  3.5°C/100m  in  the 
sector  affected  by  the  cold  feature  in  the  southeast  to 
5.0°C/100m  in  an  area  apparently  associated  with  the 
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shallower  ML D*s,  Fig  4.9  {a).  Thermocline  gradients  var 
in  a  near  meridional  manner  with  the  weaker  gradients  in 
west.  The  magnitude  of  the  gradient  varied  from  1.5°C/1 
to  1 . 8°C/100m,  Fig  4.10  (a). 

2.  Cruise  All  (29  Jun-4  Jul,  1983) 

A  general  warming  trend  was  apparent  between  the 
of  cruise  AI  and  the  beginning  of  cruise  All  with  SS 
generally  2.0°C  warmer  throughout  the  domain.  There  was  a 
a  distinct  change  in  the  SST  pattern.  Fig  4.7  (b) .  The  c 
feature  initially  observed  in  the  southeast  sector 
relaxed  southwards  and  a  second  cool  feature  was  observed 
extend  northward  into  the  southwest  of  the  domain, 
warmest  water  remained  in  the  far  northwest  sector,  wit 
near  meridional  structure  throughout  the  remainder  of 
area.  Temperatures  ranged  from  14.4°C  to  16.2°C.  The  warm 
trend  was  to  be  expected  for  the  period  involv 
Climatology  [Podeszoa,  1976]  shows  the  SST  increasing  f 
14. 1  °C  in  June  to  16.5<>C  in  July.  The  mainly  light  to  fr 
winds  and  clear  skies  prevailing  through  the  latter  part 
AI  and  early  days  of  All  could  be  expected  to  produce 
warming  trend  noted.  Satellite  imagery  once  more  agr 
well,  gualititavely,  with  the  position  of  the  major  surf 
features  shown  by  the  analysis.  Fig  3.7. 

Associated  with  the  warming  trend  and  genera 
lighter  winds,  the  BLD  shoaled  over  the  period  between 
and  All.  Layer  depths,  some  20m  shallower  than  than  th 
observed  during  AI,  were  measured  with  depths  varying  f 
near  surface  in  the  northwest  to  40m  in  the  southeast 
with  a  complex  structure  throughout  the  area. 

The  strengthening  of  the  BLG,  Fig  4.9  (b)  compa 
to  those  observed  during  cruise  AI  indicated  that 
warming  trend  had  been  limited  to  the  upper  layers  of 
ocean  and  was,  therefore,  probably  atmospherically  drii 


Magnitudes  ranged  frcm  10.0°C/100m  in  the  west  in  associa¬ 
tion  with  the  lower  SSI’s  to  3.5°C/1Q0m  in  the  centre  of  the 
domain,  loosely  associated  with  the  deeper  MLD's.  An 
increase  in  THG’s  to  tetween  1.7<>C/100m  and  2.1°C/100m  with 
a  generally  similar  pattern  to  that  of  cruise  AI,  emphasized 
the  supposition  that  the  warming  had  been  limited  to  the 
upper  layers  of  the  ocean. 

3.  Cruise  AIII  (13  Jul^iS  Jul,  1983) 

The  SST  analysis.  Fig.  4.7  (c)  showed  the  presence 
of  cool  water  along  the  southern  edge  of  the  domain  with  the 
cool  feature  in  the  southwest  during  cruise  All,  having 
relaxed  southward  and  a  further  intrusion  of  cocl  water 
extending  northward  into  the  area  from  the  southeast. 
Satellite  imagery  agreed  well  with  the  position  of  these 
features.  Fig.  3.8  .  However,  coverage  was  too  intermittent 
to  allow  the  feature  in  the  southeast  to  be  positively 
related  to  the  cool  feature  measured  in  a  similar  position 
during  cruise  AI.  Temperatures  ranged  from  12.4  to  16.0°C 
with  the  warmest  water  remaining  as  an  almost  stationary 
feature  in  the  far  northwest  of  the  domain.  A  relatively 
intense  horizontal  temperature  gradient  of  10.0°C/100m  was 
measured  in  a  30km-wide  band  around  the  cool  feature  in  the 
southeast. 

MID’s,  Fig.  4.8  (c)  were,  on  average,  deeper  than 
those  measured  during  cruise  All,  indicative  of  the  higher 
winds  prevailing  during  the  early  part  of  the  cruise 
producing  significant  mixing  in  the  surface  layers.  However, 
shoaling  of  the  MLD,  in  both  the  north  and  south  of  the 
area,  tc  10m  or  less  would  tend  to  indicate  that  the  wind 
mixing  was  restricted  to  the  area  being  sampled  by  ACANIA 
during  that  period.  The  general  horizontal  structure  of  the 
MLD  remained  complex. 


The  magnitude  of  neither  the  BLG  nor  the  THG  changed 
ouch  froa  those  measured  during  cruise  All.  However,  the 
distribution  of  the  raximum  and  minimum  BLG  was  signifi¬ 
cantly  different  from  that  of  All,  although  there  was  a  weak 
increase  in  BLG  in  association  with  the  cool  feature  in  the 
southeast  of  the  domain. 


V.  COEPABISOH  Of  OBf f BVED  ABO  MODELED  TEMPERATORB  STROCTOBB 
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For  each  of  the  parameters  discussed  in  the  Chapter  IV, 
namely,  SST,  MLD,  BLG  and  THG,  a  difference  field  was 
produced  by  subtracting  the  magnitude  of  the  parameter 
derived  from  the  model  analysis  from  that  derived  from  the 
OPTOMA5  data.  The  procedure  was  carried  out  for  each  model 
analysis  at  each  of  the  49  grid  points  in  the  domain.  The 
difference  fields  thus  calculated  were  contoured.  Figs.  5.1 
to  5.24. 

Statistics  were  calculated  from  the  contoured  difference  " 
fields  by  examining  the  percentage  of  the  total  domain 
enclosed  by  discrete  contour  intervals  in  each  difference 
field,  and  for  each  parameter.  The  ’bin-size*  for  each  field 
was  established  by  examining  the  individual  plots  and 
subjectively  selecting  an  appropriate  difference  interval 
such  that  each  bin  would  represent  between  10  and  20*  of  the 
dynamic  range  of  the  variable,  whilst  Keeping  the  area 
covered  by  each  bin  large  enough  to  make  measurement  prac¬ 
tical.  The  bin-sizes  selected  were: 

SST  -  0.4°C 
MLD  -  10.0m 
BLG  -  1 ,0°C/100m 
THG  -  0.  1oc/100m 

The  areas  within  each  bin  were  measured  using  a  Planimeter. 
Figs.  5.25  to  5.48  show  the  histograms  of  the  distributions 
obtained  by  the  above  method  and  Tables  II  to  IV  list  the 
relevant  statistics  for  each  histogram. 

A  negative  value  in  the  difference  calculations  implies 
that  the  model  analysis  was  overestimating  that  particular 
parameter  (i.e..  Difference  =  Sea  truth  -  Model  value). 
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Because  of  a  corruption  in  the  data  supplied  by  FNOC  it 
was  not  possible  to  derive  contour  fields  of  the  model  anal¬ 
yses  for  the  termination  of  cruise  AIII. 

The  larger  grid  spacing  used  in  the  TEOTS,  hemispherical 
analyses  (ca.  320ka)  meant  that  TEOTS  comparisons  were, 
necessarily,  general  in  nature  and  should  be  viewed  in  terms 
of  average  values  ,  since  it  is  unlikely  that  such  a  large 
grid  spacing  will  allow  the  analysis  to  retrieve  the  mesos- 
cale  variability  present  in  the  area  (  horizontal  vari¬ 
ability  on  a  scale  of  the  order  of  tens  of  kilometers). 

The  following  sections  discuss  the  statistics  obtained 
and  qualitatively  estimate  the  effectiveness  of  the  model 
analysis  in  representing  the  sea-truth.  A  summary  of  the 
model  performances  is  given  in  Section  A. 6  of  this  chapter. 
Throughout  the  Chapter,  sea  truth  refers  to  OPTOMA5  derived 
data* 

A.  CCHPAHISONS 
1 .  Cruise  &l 

The  analysis  dates  for  cruise  AI  were  16  June,  1983 
(initial)  and  23  June,  1983  (final)  .  FNOC  data  files  showed 
inputs  into  the  the  analyses  for  the  study  domain  of  three 
BT's  for  the  initial,  and  twenty-three  ET's  (from  R/V 
ACANIA)  for  the  final  analyses.  The  statistics  for  the 
difference  fields  derived  for  this  cruise  are  given  in  Table 
II  and  the  histograms  from  which  these  statistics  were 
derived  are  shown  in  Figs.  5.25  to  5.32 

a.  EOTS  Analysis 

The  initial  SST  analysis.  Fig.  5.1  (b)  showed  a 
nearly  meridional  structure  with  the  warmest  water  in  the 
west  of  the  domain  and  a  weakly  strengthening  horizontal 
temperature  gradient  in  the  east.  The  difference  field.  Fig. 


5.1  (c)  showed  the  analysis  to  be  overestimating  the  meas¬ 
ured  SST  for  the  majority  of  the  area.  The  mean  difference 
was  -0.8°C  and  the  maximum  difference  was  -2.0® C  associated 
with  the  intrusion  cf  cool  water  into  the  southeastern 
sector.  Fig.  5.1  (a).  The  final  analysis.  Fig.  5.1  (d) 
showed  a  general  movement  of  the  analysis  field  towards  the 
sea  truth  in  both  shape  and  magnitude.  The  shape  field  was 
only  marginaly  improved  and  continued  to  neglect  the  pres¬ 
ence  of  the  cool  feature  in  the  southeast.  The  mean  differ¬ 
ence  reduced  to  0.1 ®C  and  the  spread  of  the  difference  field 
about  this  mean  was  also  reduced  compared  to  that  of  the 
initial  analysis.  The  maximum  difference  remained  in  the 
vicinity  of  the  cool  feature  but  with  a  reduced  magnitude  of 
-0.8  °C. 

The  HLD  was  overestimated  in  both  initial  and 
final  analyses.  Fig  5.2  (b)  and  (d)  .  The  slight  shoaling  of 
the  layer  between  the  two  analyses  (2m)  had  no  significant 
effect  on  the  difference  fields.  Fig.  5.2  (c)  and  (e)  .  The 
modeled  HID  was,  on  average,  18  to  19.5m  greater  than  the 
sea  truth  and  the  maximum  difference  was  between  30  and  40m. 
The  near  uniform  structure  of  the  modeled  HLD  in  both  anal¬ 
yses  was  a  gross  misrepresentation  of  the  complex  structure 
evident  in  the  sea  truth  contour  field.  Fig.  5.2  (a). 

The  average  magnitude  of  both  the  BLG,  Fig.  5.3 
and  THG,  Fig  5.4  model  analysis  fields  also  moved  closer  to 
the  sea  truth  between  the  initial  and  final  model  analyses. 
The  BLG  analyses,  however,  showed  a  marked  difference 
between  modeled  and  sea  truth  shape  fields  for  both 
analyses. 

t.  TEOTS  Analysis 

The  initial  TEOTS  SST  analysis.  Fig.  5.5  (b) 
overestimated  the  sea  truth  in  a  similar  manner  to  that 
observed  in  the  eguivalent  EOTS  analysis  and  by 
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approximately  the  same  average  aaount  i.e. ,  a  mean 
difference  of  -0.9oc.  The  mean  difference  vas  almost  halved 
in  the  final  analysis.  Fig.  5.5  (c)  .  However,  the  error  in 
the  modeled  shape  field  remained  constant  as,  therefore,  did 
the  difference  distribution  about  the  aean.  The  difference 
fields  are  shown  in  Fig.  5.5  (d)  and  (e) . 

The  modeled  HID  was  unresponsive  to  input  data 
and  showed  a  lack  of  features.  The  TEOTS  analyses  overesti¬ 
mated  the  HID  by  an  average  of  12m.  TEOTS  HLD's  were  closer 
to  the  mean  sea  truth  HLD  than  were  the  EOTS  analyses  but, 
because  of  the  lack  cf  shape  agreement  between  model  and  sea 
truth  the  distribution  of  the  difference  field  abcut  the 
aean  was  similar  to  that  observed  in  the  EOTS  analyses. 

The  BLG,  Fig.  5.7  and  THG,  Fig.  5.8  appear  to 
respond  to  the  input  data  by  moving  closer  to  the  sea  truth 
in  the  mean.  However,  the  lack  of  shape  agreement  between 
model  and  sea  truth  contour  fields  remained  unchanged.  The 
TEOTS  analysis  was  initially  further  from  the  sea  truth  and 
responded  more  slowly  to  input  data  than  the  EOTS  analysis. 

2 .  Cruise  £11 

The  analysis  dates  for  cruise  All  were  30  June,  1983 
(initial)  and  7  July,  1983  (final).  The  inputs  to  the  anal¬ 
yses  were  3  BT's  and  nil  for  the  initial  and  final  analyses, 
respectively.  Statistics  for  the  model/sea  truth  comparisons 
are  given  in  Table  III  and  the  histograms  from  which  the 
statistics  were  derived  are  shown  in  Figs.  5.33  to  5.40 

a.  EOTS  Analysis 

The  initial  SST  analysis.  Fig.  5.9  (b)  showed 
similar  shape  characteristics  to  those  observed  in  the  23 
June  EOTS  analysis  (cruise  AI)  and  a  general  warming  trend 
of  0.5OC  was  indicated.  Over  the  same  period  the  sea  truth 
analysis.  Fig.  5.9  (a)  indicated  an  increase  in  SST  of  2.0°C 


•1  ya. 


and  a  marked  change  in  the  shape  of  the  contour  field.  The 
initial  difference  field.  Fig.  5.9  (c)  showed  the  analysis 
to  be  an  underestimate  of  the  sea  truth  by  an  average  of 
1.0°C  with  a  maximum  difference  of  1.6«c  and  a  general  lack 
of  agreeient  between  nodel  and  sea  truth  shape  fields  was 
apparent. 

The  final  model  analysis.  Fig.  5.9  (d) ,  showed  a 
trend  towards  a  meridional  SST  structure  (away  from  the 
zonal  structure  apparent  in  the  sea  truth  field)  with  no 
change  in  the  average  SST.  The  mean  difference  remained 
similar  to  that  calculated  for  the  initial  analysis. 
However,  because  of  the  increased  disagreement  between  the 
model  and  sea  truth  shape  fields,  the  spread  of  the  differ¬ 
ence  distribution  about  the  mean  increased. 

The  FNOC  climatological  depth/teaperature 
profile  for  the  June/July  period  in  the  region  of  the  study 
domain  predicts  a  SST  of  14.2°C  which  was  in  close  agreement 
with  the  mean  SST  of  14.4°C  observed  in  the  model  analysis. 

The  initial  HLO  analysis.  Fig.  5.10  (b)  showed 
no  significant  change,  in  shape  or  magnitude  from  that 
observed  in  the  23  June  analysis.  Fig.  5.2  (d).  However, 
changes  in  the  sea  truth  field.  Fig.  5. 10  (a)  ,  between 
cruises  AI  and  All  meant  that  the  difference  field.  Fig. 
5.10  (c)  ,  was  markedly  different  from  that  of  the  23  June 
analysis.  Although  the  final  HLD  analysis.  Fig  5.10  (d)  did 
indicate  a  slight  shoaling  (5s  to  10m) ,  the  HLD  was  overes¬ 
timated  by  the  model  analyses  throughout. 

The  shoaling  of  the  HLD  between  the  initial  and 
final  analyses  noted  above  is  compatible  with  the  model 
analysis  moving  towards  the  FNOC  July  climatology  for  the 
region  which  shows  a  change  in  HLD  from  55m  in  June  to  37m 
in  July. 

The  model  BLG's  were  overestimated  and  the  shape 
of  the  contour  field  in  poor  agreement  with  that  of  the  sea 


truth.  A  smoothing  of  the  modeled  BLG  field  to  a  near 
constant  value  in  the  final  analysis  did,  however,  reduce 
the  mean  difference  from  the  1.2°C/100m  observed  in  the 
initial  analysis  to  0.3°C/100m.  The  results  of  the  compari¬ 
sons  are  shown  in  Fig.  5.11. 

Model  and  sea  truth  THG's  were  in  good  agreement 
for  both  analyses.  Fig.  5.12. 

b.  TEOTS  Analysis 

The  warming  trend  noted  in  the  EOTS  analyses  was 
also  observed  in  the  TEOTS  analyses.  Fig.  5. 13  The  mean 
difference  decreased  from  0.40c  in  the  initial  TEOTS  anal¬ 
ysis  tc  near  zero  in  the  final  analysis.  The  horizontal 
structure  of  the  TEOTS  analyses  remained  meridional  with  a 
weakly  strengthening  horizontal  temperature  gradient  in  the 
east  of  the  domain  compared  to  the  complex  near  zonal  struc¬ 
ture  observed  in  the  sea  truth  SST. 

Modeled  MID*s,  Fig.  5.14  remained  unchanged  by 
any  significant  amount  in  both  magnitude  and  structure  from 
those  of  the  TEOTS  analyses  of  cruise  AI.  The  modeled  MLD* s 
were  overestimated  by  between  5m  and  35m,  when  compared  to 
the  sea  truth. 

The  BLG,  Fig.  5.15  was  overestimated  by  the 
model  in  both  the  initial  and  final  analyses  with  the 
magnitude  of  the  mean  difference  increasing  by  the  final 
analysis. 

The  THG,  Fig.  5. 16  was  in  reasonably  good  agree¬ 
ment  with  the  sea  truth  throughout. 

3 •  Cruise  A  III 

As  previously  stated,  model  analyses  were  not  avail¬ 
able  for  the  termination  of  this  cruise.  There  were  no  BT 
inputs  to  either  model  for  the  10  July,  1983  analysis. 


The  only  significant  change  in  either  model  was  a 
general  increase  in  the  SST  {0. 5°C)  indicated  in  the  ECTS 
analysis. 

The  significant  changes  in  the  difference  fields. 
Figs.  5.17  to  5.24  were  due  to  changes  in  the  sea  truth 
fields  between  cruises  All  and  AIII. 

The  statistics  of  the  difference  fields  and  the 
relevant  histograms  are  shown  in  Table  IV  and  Figs.  5.41  to 
5.48  As  with  cruise  All,  the  TEOTS  SST  and  MLD  analyses 
were,  on  average,  in  closer  agreement  with  the  sea  truth 
than  the  EOTS  analyses. 

4.  Effectiveness  of  the  EOTS  Analysis  in  Representing 

Input  Data 

To  estimate  the  accuracy  of  the  EOTS  analysis  in 
responding  to  input  data  the  author's  analysis  program  was 
revised  to  operate  on  a  4x4  grid  with  a  40km  grid  spacing, 
similar  to  that  used  in  the  EOTS  analysis  program.  The 
twenty-three  BT's  input  to  the  EOTS  23  June  analysis  were 
used  as  input.  The  resulting  SST  contour  field  and  that 
produced  by  the  equivalent  EOTS  analysis  are  shown  in  Fig. 
5.49. 

From  Fig.  5.49  it  can  be  seen  that  the  EOTS  analysis 
did  net  respond  accurately  to  the  input  data.  The  accuracy 
of  the  author's  analysis  program  had  previously  been  veri¬ 
fied  by  comparisons  with  manually  contoured  SST  fields  with 
favourable  results.  The  reason  for  the  lack  of  resolution  in 
the  ECTS  analysis  is  not  fully  understood  but  is  probably 
due  to  excess  weighting  on  the  first  guess  field  used  in  the 
EOTS  FIB  analysis.  This  first  guess  field  provides  the  shape 
information  to  which  any  input  data  is  fitted  in  a  least 
squares  sense. 


5*  Validity  Qu^si-Synoptic  Assumption 

lack  of  synopticity  in  the  0PT0NA5  data  proved  to  be 
unimportant  since  the  only  inputs,  of  any  significance,  to 
the  F HOC  models  were  the  twenty-seven  BT  submissions  prior 
to  the  23  June  analyses.  The  fact  that  the  BT's  were  input 
to  the  models  as  one  data  set  made  them  synoptic  to  the 
computer  model  and,  therefore,  a  synoptic  assumption  is 
valid  for  the  sea  truth  analysis  with  which  the  model  anal¬ 
yses  were  compared. 

6 .  Summary 

a.  SST 

The  EOTS  SST  analysis  responded  to  the  magnitude 
changes  dictated  by  input  data  (when  sufficient  data  was 
available)  but  horizontal  structure  was  not  resolved  on  a 
scale  commensurate  with  the  resolution  capability  of  the 
40km  grid  spacing  used.  In  the  absence  of  input  data  the 
EOTS  SST  moved  towards  climatology  and  became  even  less 
structured. 

There  was  a  measureable  response  to  input  data 
in  the  TEOTS  SST  analysis.  The  magnitude  of  the  response  was 
less  than  that  observed  in  the  EOTS  analysis  and  the  lack  of 
agreement  in  the  shape  field  between  model  and  sea  truth  was 
even  more  pronounced,  as  might  have  been  expected  frcm  the 
larger  grid  spacing  used  in  the  TEOTS  analysis.  TECTS 
perfomed  better,  in  a  general  sense,  than  EOTS  when  data 
were  not  available  and  the  move  towards  climatology  noted  in 
the  ECTS  analysis  was  not  observed. 

b.  NLD 

Neither  aodel  produced  a  realistic  analysis  of 
the  sea  truth  NLD  whether  input  data  were  available  or  not. 
Both  analyses  tended  to  overestimate  the  magnitude  of  the 


HID,  with  the  TEOTS  analysis  giving  a  consistently  better 
estimate  of  the  sea  truth  HLD  than  the  EOTS  analysis. 

c.  BLG  and  THG 

The  lack  of  agreenent  between  modeled  and  sea 
truth  HID  made  comparisons  of  BLG  and  THG  difficult  (Both 
parameters  were  defined  in  terms  of  the  layer  depth).  Any 
agreement  noted  between  these  two  parameters  was  therefore 
regarded  as  purely  coincidental. 

B.  INTEB-ARIUAL  VARIABILITY 

The  mesoscale  variability  of  the  study  domain  has  been 
dealt  with  in  some  detail  in  the  previous  sections. 
Underlying  this  mesoscale  variability  is  a  variability  with 
a  time  scale  of  the  order  of  years  (the  effect  of  El  Nino) . 
The  region  was  known  to  be  influenced  by  such  an  event 
during  the  period  of  the  present  study.  To  illustrate  the 
effect  of  such  an  anomaly,  mean  depth/temperature  profiles 
were  compared  for  July,  1983  and  July/August,  1982  in  the 
study  domain.  The  mean  profiles  were  also  compared  to  FNOC 
July  climatology  for  the  region.  The  mean  profile  for  1982 
was  derived  from  XE1  data  collected  during  the  0PT0HA2 
cruise  of  July/August,  1982  (The  area  of  the  CCS  sampled  by 
this  cruise  was  very  close  to  that  sampled  in  the  0PT0MA5 
series  of  cruises)  and  that  for  1983  was  derived  from  XBT 
data  collected  during  cruise  AIII  of  the  0PT0HA5  series 
(July,  1983).  Figs.  5.50  to  5.52  show  the  results  of  the 
comparisons. 

FNOC  climatology  was  in  reasonable  agreement  with  the 
1982  mean  profile  (assumed  non-El  Nino  year).  Fig  5.50  .  Ihe 
major  discrepancy  between  climatological  and  measured 
profiles  occurred  in  the  upper  part  of  the  water  column 
(surface  to  50m)  which  could  be  expected  since  this  is  the 


region  subject  to  maximum  atmospheric  forcing  and,  there¬ 
fore,  the  short  time  scale  variability  impossible  to  predict 
with  climatology. 

The  presence  of  a  thermal  anomaly  in  the  depth/ 
temperature  profile  was  clearly  evident  when  the  1982 
profile  was  compared  to  that  from  1983,  Fig.  5.51.  The 
assumed  £1  Nino  event  was  present  as  a  layer  of  warmer 
waters  between  the  near  surface  and  approximately  200m  with 
a  maximum  anomaly  of  2. 5°C  at  a  depth  of  60m.  The  inability 
of  a  static  climatology  to  represent  such  an  event  was 
demonstrated  by  comparing  the  1983  mean  depth/temperature 
profile  with  that  obtained  from  FNOC  climatology.  Fig.  5.52. 

To  obtain  a  measure  of  the  operational  significance  of 
an  anomalous  event  such  as  the  one  described  above  transmis¬ 
sion  loss  curves  were  computed  from  the  FACT9H  transmission 
loss  model  for  theJuly  1983  climatological  and  mean  depth/ 
temperature  profiles  .  The  only  measurable  effect.  Figs. 
5.53  to  5.55  was  a  reduction  in  convergence  zone  range  in 
the  curves  derived  ficm  the  climatological  profiles  (approx¬ 
imately  3%)  and  a  decrease  in  convergence  zone  gain 
(approximately  5db)  . 

C.  RIS01TS  OF  RECENT  FNOC  STUDIES 

Recent  modif ications  have  been  made  to  both  EOTS  and 
TEOTS  analysis  models  (Frost,  pers.  comm.)  and,  in  partic¬ 
ular,  the  model  sensitivity  to  input  changes  in  PLD  (MLD) 
has  been  modified.  The  incorporation  of  the  TEOTS  analysis 
into  the  finer  grid  regional  model  has  also  been  achieved 
for  certain  areas.  The  effect  of  these  modif ications  is  the 
subject  of  an  ongoing  investigation  by  FNOC,  based  on  inputs 
from  the  Bay  of  Biscay  and  the  Western  Mediterranean  Sea 
regional  analyses. 


The  initial  results  from  the  study  are  promising  with 
TEOTS  apparently  showing  considerable  skill  in  predicting 
changes  in  layer  depth  even  in  the  absence  of  input  data 
(Frost,  pers.  comm.).  No  information  is  available  from  this 
study  related  to  the  skill  of  the  analysis  in  resolving  the 
horizontal  structure  in  the  region. 


a.  Cruise  AI  Analysis 


b.  EOTS  Analysis  16  Jun 


c.  Difference  Field  16  Jun 


Figure  5.2  Contoured  BID  Fields  for  EOTS/Cruise  AI 

(contour  units  b)  . 
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a.  Cruise  AI  Analysis 


Figure  5.4  Contoured  THG  Fields  for  EOTS/Cruise  AI 
(contour  units  °C/1000n). 
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a.  Cruise  All  Analysis 


Figure  5.9  Contoured  SST  Fields  for  EOTS/Cruise  ill 

(contour  units  °Cx10) . 
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a.  Cruise  All  Analysis 


Figure  5.11  Contoured  BLG  Fields  for  EOTS/Cruise  All 
(contour  units  °C/1000n)  . 
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b.  TEOTS  Analysis  30  Jun 


c.  Difference  Field  30  Jun 


a.  Cruise  All  Analysis 
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d.  TEOTS  Analysis  7  Jul 


e.  Difference  Field  7  Jul 
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Figure  5.15  Contoured  BLG  Fields  for  TEOTS/Cruise  A 

(contour  units  °C/1000) . 
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d.  TEOTS  Analysis  7  Jul 


Figure  5.16  Contour* 
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a.  Cruise  AIII  Analysis 


Pigure  5.18  Contoured  HLD  Fields  for  EOTS/Cruxse  AIII 

(contour  units  ■) . 
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a.  Cruise  AIII  Analysis 


Figure  5.19  Contoured  BLG  Fields  for  EOTS/Cruise  AIII 

(contour  units  °C/1000) . 
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a.  Cruise  AIII  Analysis 


Figure  5.20  Contoured  TBG  Fields  for  EOTS/Cruise  AIII 

(contour  units  oc/1000) • 
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a.  Cruise  AIII  Analysis 


Pigure  5.21  Contoured  SST  Fields  for  TEOTS/Cruise  AIII 

(contour  units  °Cx10) . 

87 


IN Ji'Wli J/109  IV  03Dna04id3iV 


«lt 


VI ' 


a.  Cruise  AIII  Analysis 


Figure  5.23  Contoured  BLG  Fields  for  TEOTS/Cruise  AIII 

(contour  units  °C/1000) . 
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a.  Cruise  AIII  Analysis 


Figure  5.24  Contoured  THG  Fields  for  TBOTS/Cruise  Alii 

(contour  units  °C/1000) . 
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Figure  5.43  BL6  Difference  Distribution  EOTS/Cruise  ill 

10  Jul  Analysis. 
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5.53  Transmission  loss  Curves-S/S  Disposition 


Climatology  b.  0PT0MA5  AIII  (mean  profile) 
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Figure  5.54  Transmission  loss  Curres-S/D  Disposition 
a.  Cliaatology  b.  OPTOH15  AIII  (aean  profile) . 
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TABLE  II 

CoaparisoD  Statistics  Craise  AI 


i 

PARAMETER 

i  . 

MODEL 

EOTS 

TEOTS 

ANALYSIS  DATE 

16  Jun 

23  Jun 

16  Jun 

23  Jun 

SST 

mean 

-0.80 

0.10 

-0.90 

-0.50 

std.  deviation 

0.50 

0.30 

0.50 

0.40 

(•  C) 

med ian 

-0.60 

0.20 

-1.00 

-0.60 

skewness 

-0.40 

-0.30 

-0.30 

-0.20 

MLD 

mean 

warn 

-17.9 

-6.50 

-5.70 

std.  deviation 

7.70 

7.30 

7.70 

7.70 

(m) 

med ian 

-15.0 

-15.0 

-5.00 

-5.00 

skewness 

0.10 

0.20 

-0.10 

-0.20 

! 

BLG 

mean 

-1.26 

-0.05 

-3.24 

-1.50 

std.  deviation 

1.16 

0.70 

0.90 

1.00 

(•C/lOOm) 

aed ian 

-0.50 

-0.50 

-3.50 

-1.50 

skewness 

-0.80 

0.50 

-0.10 

-0.40  : 

| 

TKG 

mean 

L - - 

0.08 

0.02 

-0.11 

- 1 

-0.09 

std.  deviation 

mgm 

0.07 

0.12 

1 

0.09 

(°C/100m) 

med ian 

0.05 

0.05 

-0.05 

-0.05  j 

skewness 


0.50 


0.30 


-0.60 


-0.60 


TABLE  III 

Coaparison  Statistics  Cruise  All 


MODEL 

EOTS 

TEOTS 

PARAMETER 

ANALYSIS  DATE 

30  Jun 

07  Jul 

30  Jun 

07  Jul 

mean 

0.90 

1.03 

0.40 

-0.04 

sCd.  deviation 

0.40 

0.50 

0.60 

0.60 

SST 
(°  C) 

median 

1.00 

1.00 

0.60 

0.20 

skewness 

-0.90 

-0.70 

-0.02 

0.01 

mean 

-30.9 

-26.6 

-11.2 

-17.4 

MLD 

std.  deviation 

8,30 

7.70 

16.8 

7,89 

(m) 

median 

-25.0 

-25.0 

-15.0 

-15.0 

skewness 

-0.20 

-0.60 

2.16 

-0.80 

mean 

-1.20 

-0.30 

-2.80 

-4.24 

BLG 

(•C/lOOm) 

std.  deviation 

1,90 

1.40 

1.60 

1.49 

med ian 

-1.50 

-0.50 

-3.50 

I 

O 

skewness 

0.90 

1.40 

1.20 

0.90 

mean 

0.09 

0.25 

0.08 

-0.04 

THG 

("C/lOOm) 

std.  deviation 

0.07 

0.12 

0.11 

0.11 

med ian 

0.20 

0.30 

0.15 

-0,05 

skewness 

-0.60 

-1.00 

-0.30 

-0.20 

TABLE  IT 

Coaparieca  Statistics  Cruise  AIII 


parameter 

MODEL 

EOTS 

TEOTS 

ANALYSIS  DATE 

10  Jul 

21  Jul 

10  Jul 

21  Jul 

mean 

0.50 

0.12 

SST 

std.  deviation 

0.60 

0.60 

(•  C) 

med ian 

0.60 

0.20 

skewness 

0.20 

0.10 

mean 

-5.60 

1.50 

MLD 

std.  deviation 

16.90 

13.10 

(m) 

median 

-5.00 

5.00 

skewness 

-0.40 

1 

O 

o 

mean 

0.05 

-3.90 

BLG 

C*C/100m) 

std.  deviation 

1.33 

1.15 

median 

0.50 

-3.50 

skewness 

0.20 

-0.05 

mean 

0.20 

0.04 

THG 

std.  deviation 

0.14 

0.14 

(•C/lOOm) 

med ian 

0.20 

0.05 

skewness 

-0.20 

0.30 

t 
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▼I.  iCOOSTIC  IBP ACT  21  ESHfiES  IN  FNOC  OCEAN  THEBHAL 

H&OC12A1  AMALjgis 


A.  DEFINITION  OF  ACOOSTIC  TEENS 

•  Convergence  Zone  (CZ)  -  Region  of  focussing  or  convergence 
cf  sound. 

•  Direct  Range  -  The  shortest  distance  from  the  receiver  at 
which  the  transmission  loss  exceeds  the  Figure  of  Merit. 

•  Figure  of  Merit  (FOM)  -  The  aaxiaun  allowable  one  way 
transaission  loss  (passive  sonars)  for  which  a  50%  prob¬ 
ability  of  detection  is  possible  for  a  particular  sonar. 

•  Leakage  -  Loss  of  energy  from  a  surface  duct  to  the  main 
sound  channel. 

•  Low-Frequency  Cut-Off  -  Frequency  below  which  no  trapped 
normal  mode  may  exist  for  a  given  channel  or  duct. 

•  Sound-Speed  Profile  (SSP)  -  Depth/sound-speed  combinations 
defining  the  vertical  variation  in  sound  speed. 

•  Transmission  Loss  (TL)  -  Acoustic  parameter  which  quanti¬ 
tatively  describes  the  weakening  of  sound  intensity  between 
a  point  1m  from  the  source  and  a  point  at  a  distance  in  the 
ocean. 

•  Vertical  Beam-Width  -  Included  angle  between  the  maximum 
upward  and  downward  propagating  rays. 

B.  ASSUMPTIONS 

The  cnly  variable  in  the  range  calculations  was  to  be 
the  SSP.  To  achieve  this  goal  certain  representative, 
acoustic  parameters  were  assumed  constant  throughout  and 
were  used  as  standard  inputs  to  all  transmission  loss 
calculations. 

•  The  ocean  bottom  was  assumed  flat. 


•  A  standard  bottom- less/grazing  angle  curve  applied  to  all 
frequencies. 

•  A  flat  sea  surface  (sea  state  1)  input  to  the  FACT  models 
(assused  flat  in  the  FE  aodel) 

•  The  FOH*s  used  in  the  range  calculations  were  80db  and 
90db.  These  figures  are  representative  of  the  range  of  FOB’S 
encountered  in  operational  systems. 

•  Source/receiver  dispositions  used  were: 

20a/20a-shallow/shallow  (S/S) 
20a/200a-shallow/deep  (S/0) 
200a/200a-deep/deep  (0/0) 

The  S/S  disposition  was  intended  to  give  a  neasure  of 
the  surface  duct  acoustic  propagation;  S/D  of  cross-layer 
propagation  and  D/D  of  deep  sound  channel  propagation.  All 
three  coabinations  can  be  utilized  by  presently  operational, 
shipberne  systems. 

C.  DERIVATION  OF  SOMD-SPBED  PROFILES 

Depth/temperature  profiles  were  constructed  for  both  the 
EOTS  and  OPTOHA5  analyses  at  the  positions  of  aaxinua  error, 
between  the  two  analyses,  for  SST,  NLD,  BIG  and  THG.  The 
position  of  the  maximum  error  was  obtained  from  the  differ¬ 
ence  field  discussed  in  Chapter  V  and  the  depth/temperature 
profiles  were  generated  by  utilizing  the  contour  fields 
derived  as  described  in  Chapter  IV. 

SSP's  were  computed  from  the  temperature  profiles  using 
Rilson’s  Equation.  The  depth/salinity  information  was 
obtained  from  the  DE  STEIGEOR,  deep  CTO  data  (Chapter  II) 

Only  the  data  relevant  to  cruises  AI  and  All  were 
utilized.  Long  range  prediction  (CZ)  comparisons  were  made 
between  profiles  based  on  the  0PT0HA5,  EOTS  and  TEOTS  anal¬ 
yses.  Direct  range  comparisons  were  limited  to  profiles 
based  on  the  0PT0HA5  and  EOTS  analyses  only. 


126 


The  SSP's  derived  for  each  aodel  analysis  (and  all 
subsequent  acoustic  calculations  based  on  those  SSP's)  were 
designated  according  to  the  paraneter,  and  analysis  date  eg. 
SST1  refers  to  a  profiles  derived  at  the  position  of  aaxiaua 
difference  between  acdeled  and  0PT0HA5  SST  for  the  initial 
aodel  analysis  of  the  particular  0PT0HA5  cruise  being 
studied;  MLD2  refers  to  the  profile  derived  at  the  position 
of  aaxinua  difference  in  HLD's  for  the  final  aodel  analysis 
of  the  cruise,  etc. 

D.  TBAHSHISSIOI  LOSS  BODEL  COMPABISOHS 

Differences  in  aodel  physics  nade  direct  comparisons 
between  FACT9H  and  the  PE  Models  unwieldy,  as  did  limita- 
tions  in  the  aodel  computer  programs.  One  of  the  major 
differences  was  that  the  PE  aodel  has  a  vertical  bean width 
limitation  which  causes  a  program  termination  if  exceeded. 
The  limitation  is  a  function  of  the  SSP  and,  for  the 
profiles  utilized  in  this  study,  the  maximun  vertical  bean- 
width  which  could  be  maintained  was  20°.  Mo  such  limitation 
exists  fer  the  FACT9H  aodel  which  suas  acoustic  intensities 
froa  all  forward  directions,  equivalent  to  a  vertical  bean- 
width  of  1800.  The  effect  of  this  beaawidth  liaitation  is 
that  the  PE  aodel  will  predict  consistently  greater  trans¬ 
mission  loss  for  a  given  range  than  will  the  FACT9H  model, 
assuming  all  other  inputs  to  each  model  are  identical. 
However,  the  effect  of  SSP  changes  relating  to  temperature 
profile  differences  between  modeled  and  0PT0HA5  ocean 
thermal  structures  were  computed  on  a  model-by-model  basis 
and  the  results  compared.  Table  7.  An  example  of  the  trans¬ 
mission  loss  curves  computed  from  the  two  models  for  the 
same  SSF,  Fig.  6.1,  clearly  shows  the  effect  of  the 
differing  model  physics. 


Table  7  shows  the  Parameter  Error  (difference  between 
OPTOHA5  and  EOTS  values  for  the  parameter,  expressed  as  a 
percentage  of  the  OPTCHA5  value);  the  Range (B)  (direct  range 
(ka)  coaputed  froa  the  transaission  loss  curve  based  on  the 
OPTO HAS  data) ;  Range  Error  (E)  (difference  between  direct 
range  obtained  utilizing  0PT0MA5  data  and  that  utilizing 
EOTS  analysis  data,  expressed  as  a  percentage  of  R).  A 
negative  sign  indicates  that  the  nodel  analysis  had  produced 
an  overestiaate  of  that  particular  paraaeter.  The  compari¬ 
sons  shown  in  Table  V  are  for  a  600Hz  frequency. 

The  lack  of  agreeaent  between  the  two  transaission  loss 
aodels  was  apparent  with  the  PE  nodel  predicting  signifi¬ 
cantly  shorter  ranges  than  the  FACT9H  nodel.  The  PE  nodel 
was  also  far  aore  sensitive  to  variations  in  the  input  SSP, 
particularly  when  both  source  and  receiver  were  shallow,  and 
showed  a  dependence  on  the  SSP  for  the  S/D  and  D/D  disposi¬ 
tions,  which,  as  discussed  in  later  sections,  was  not  the 
case  for  the  FACT9H  aodel.  Convergence  zone  ranges  coaputed 
by  both  FACT9H  and  the  PE  aodel  were  in  good  agreement. 
However,  the  gain  in  acoustic  energy  at  the  convergence 
zones  was  approximately  5db  lower  in  the  PE  model  computa¬ 
tions  than  those  of  FACT9H. 

The  remaining  comparisons  discussed  in  this  chapter  were 
based  on  the  FACT9H  transaission  loss  computations. 

E.  SUBBART  OF  ACOUSTIC  COHPUTATIOHS 

When  comparing  ranges  predicted  utilizing  measured  and 
modeled  profiles,  an  error  of  less  than  20%  was  assumed 
operationally  insignificant.  This  assumption  was  believed 
reasonable  when  the  uncertainties  associated  with  the  FACT9H 
acoustic  aodel  accuracy  (eg.,  the  lack  of  attenuation  due  to 
absorption  and  volumetric  scattering)  are  considered. 


1.  Ism  Raffle  (£Z)  E£§Al£U2&s 

Comparisons  of  CZ  ranges  computed  using  SSP’s 
derived  from  OPTOMA5,  EOTS  and  TEOTS  data.  Tables  VI  to  IX 
revealed  a  relative  lack  of  sensitivity  of  the  CZ  range  to 
thermal  structure  variations  of  the  order  of  magnitude 
present  between  modeled  and  sea  truth  analyses.  The  average 
CZ  range  computed  from  sea  truth  data  was  50 km,  varying 
between  44km  and  62km.  Ranges  computed  utilizing  SSP's 
derived  from  FNOC  model  analyses  were  within  10*  cf  those 
computed  utilizing  sea  truth  data.  The  maximum  errors  were 
observed  when  both  source  and  receiver  were  shallow. 

2.  Direct  Range  Comparisons 

Subjective  analysis  of  the  transmission  loss  curves 
computed  for  the  long  range  comparisons  discussed  above, 
combined  with  an  understanding  of  ray-theory  physics,  led  to 
the  direct  range  comparisons  being  concentrated  on  the 
surface  duct  propagation  paths  (i.e.,  shallow  source  and 
receiver).  A  limited  number  of  comparisons  were  made  between 
direct  ranges  obtained  from  the  transmission  loss  curves 
computed  utilizing  OPTOHA5  and  those  utilizing  EOTS  derived 
profiles  for  cross-layer  and  below  layer  propagation  paths. 
The  results  from  these  comparisons.  Tables  X  and  XI  confirm 
the  supposition  that  the  FACT9H  computed  transmission  losses 
for  the  S/D  and  D/D  propagation  paths  are  relatively  insen¬ 
sitive  to  SSP  variations  of  the  order  of  magnitude  present 
between  profiles  derived  from  0PT0MA5  data  and  those  from 
the  EOTS  analyses.  In  the  remainder  of  this  discussion 
’range  error*  will  refer  to  the  difference  between  ranges 
calculated  from  transmission  loss  curves  computed  from  sea 
truth  profiles  and  those  computed  from  EOTS  profiles. 

The  maximum  range  error  for  the  cross-layer  propaga¬ 
tion  path  was  18*  on  a  range  of  8km,  with  the  majority  of 


the  range  errors  within  1 OX  of  a  sea  truth  range  of  8km  at 
the  50Hz  frequency.  However,  a  weak  frequency  dependence  was 
observed  and  errors  cf  the  order  of  15%  on  a  range  of  8kn 
were  measured  at  600Hz,  associated  with  the  lover  FOH.  The 
maximum  error  also  occurred  at  this  frequency,  associated 
with  the  lower  FOH,  Table  X  . 

The  range  errors  measured  for  the  below  layer  propa¬ 
gation  paths  were  generally  negligible  (<5%)  and  independent 
of  frequency  and  FOH,  Table  XI  .  However,  for  profiles 
SSTIand  SST2  errors  of  the  order  of  15%  on  a  range  of  18km 
for  the  lower  FOH  and  15%  on  a  range  of  77kn  for  the  higher 
FOH,  were  measured.  The  presence  of  a  weak  shallow  sound 
channel  in  the  sea  truth  SSP  is  believed  to  be  the  cause  of 
these  significantly  higher  range  errors  and  would  also 
explain  the  longer  direct  ranges  associated  with  the  sea 
truth  profile  (i.e.,  18km  compared  to  the  average  of  8km  for 
the  remaining  profiles) 

For  the  surface  duct  propagation  paths  (S/S) ,  trans¬ 
mission  loss  curves  were  computed  at  frequencies  of  50,  300 
600  and  1000Hz  for  profiles  associated  with  cruise  Al  and  at 
300,  600  and  1000Hz  for  those  associated  with  cruise  All. 
The  50Hz  frequency  was  removed  from  the  comparisons  in 
cruise  All,  Table  XIII,  once  the  range  errors  were  shown  to 
be  negligible  at  this  frequency  for  the  cruise  AI  compari¬ 
sons,  Table  XII  . 

From  the  range  comparisons.  Table  XII  and  XIII  it 
can  be  seen  that  at  the  higher  FOH  and  for  frequencies  below 
1000Hz,  the  transmission  loss  is  insensitive  to  SSP  varia¬ 
tions  of  the  order  of  magnitude  observed  with  a  range  error 
of  less  than  10%  on  a  mean  range  of  55km. 

For  range  calculations  based  on  the  lower  FCH, 
significant  range  errors  were  observed  for  frequencies  of 
300Hz  and  above.  The  range  error  was  also  observed  to  be  a 
function  of  frequency.  Plots  of  range  error  versus  parameter 


error  for  the  four  parameters  being  studied  revealed  a 
strong  dependence  of  the  range  error  on  the  MLD  error.  Fig. 
6.2  .  The  curves  shown  in  Fig.  6.2  were  fitted  to  the  points 
in  an  RMS  manner  and  do  not  imply  a  functional  relationship 
between  HID  and  range  error  but  are  intended  to  shew  a 
general  trend.  The  scatter  of  the  points  in  Fig.  6.2  indi¬ 
cate  the  complex  interaction  between  the  various  parameters. 
From  the  curves  it  can  be  seen  that  the  range  error 
increased  with  increasing  MLD  error  and  frequency,  for 
freguecies  above  approximately  300Hz.  The  presence  cf  a 
non-zero  range  error  for  a  zero  MLD  error  was  a  further 
indication  of  the  effect  of  the  other  parameters.  Piets  of 
range  error  versus  errors  in  the  remaining  three  parameters 
were  also  generated.  Figs.  6.3  to  6.5  These  curves  appeared 
to  indicate  a  marked  dependence  of  range  error  on  each 
parameter.  However,  plots  of  MLD  error  versus  the  errors  in 
the  remaining  parameters.  Figs.  6.6  to  6.8  revealed  that  the 
major  range  error  dependence  was  indeed  on  MLD  error;  e.g. 
the  trend  in  range  error  when  compared  to  SST  error  is  in 
the  same  general  sense  as  the  trend  in  MLD  error  when  that 
too  is  compared  to  SSI  error. 
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Figure  6.2  Plot  of  HLD  Error  vs.  Range  Error. 
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Figure  6.3  Plot  of  SST  Error  vs.  Hange  Error. 
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Figure  6.4  Plot  of  BIG  Error  ts.  Range  Error. 
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Figure  6.6  Plot  of  SST  Error  vs 
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Figure  6.7  Plot  of  BIG  Error  vs.  HID  Error. 
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Figure  6.8  Slot  of  THG  Error  ts.  HLD  Error. 
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TABLE  VII 

TEOTS/Croise  Ai  CZ  Range  Coaparisons 
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TABLE  IX 

TEOTS/Cruise  All  CZ  Range  Comparisons 


TABLE  X 

Cross-Layer  Propagation  Conparisons  (S/D) 


TABLE  XI 

Belov-Layer  Propagation  Comparisons  (D/D) 


THC2 


TABLE  XII 

Surface  Duct  propagation  Coaparisons  (S/S) -Cruise  AI 


•K 


d3(j 


THG2 


VII.  S0HHA8T  ASP  CONCLUSIONS 

The  grid  spacing  used  in  the  FHQC  ocean  thermal  analysis 
models  {i.e.,  40km  for  EOTS  and  320km  for  TEOTS)  made  it 
extremely  unlikely  that  they  would  be  capable  of  recovering 
the  mesoscale,  horizontal  variability  (of  the  order  of  30  to 
40km)  dominating  the  study  domain. 

Within  the  limitations  of  the  grid  spacing  used,  the 
EOTS  SST  analysis  responded  well  to  input  data.  However,  the 
analysis  methodology  is  believed  to  have  smoothed  the  input 
field  to  the  extent  that  significant  detail  was  lost. 

The  hemispherical  TEOTS  SST  analysis  responded  to  data 
input  in  magnitude  but  was  totally  insensitive  to  the  hori¬ 
zontal  shape  information  provided  by  the  input  data.  This 
lack  of  response  is  believed  to  be  a  function  of  the  coarse 
grid  spacing  of  the  TEOTS  analysis  field. 

Both  EOTS  and  TEOTS  were  insensitive  to  input  mixed 
layer  depth  (HLD)  information.  The  EOTS  HLD  analysis 
remained  almost  constant  throughout  the  study  period  in  both 
magnitude  and  shape.  The  TEOTS  HLD  analyses  did  change  over 
the  study  period  and  were,  on  average,  in  closer  agreement 
with  the  measured  HLD  magnitudes,  but  had  an  almost  constant 
shape  field.  The  lack  of  detail  in  both  EOTS  and  TEOTS  HLD 
analyses  was  in  marked  contrast  to  the  complex  and  variable 
measured  HLD  shape  fields.  The  improved  performance  of  the 
TEOTS  analysis  over  the  EOTS  analysis  is  believed  to  be  the 
result  of  the  influence  of  the  TOPS  element  on  the  analysis. 

The  response  of  the  models  to  data  input  for  depths 
below  the  HLD  was  inconclusive  since  the  parameters  tested 
(BLG  and  THG)  were  functions  of  the  HLD  and,  thus,  errors  in 
the  HLD  analysis  would  propagate  into  these  remaining 
fields. 


The  FNOC  procedure  of  using  a  single  climatological 
depth/teaperature  profile  for  areas  such  as  the  study  domain 
vas  shown  to  be  inaccurate.  The  study  domain  was  shown  to 
wary  on  a  secular  time  scale  (El  Nino  effect)  with  a  2  to 
3°C  thermal  anomaly,  compared  to  FNOC  climatology,  measured 
at  depths  of  50  to  100a  during  the  study  period.  The  effect 
of  a  poor  climatology  would  be  reflected  in  the  model  anal¬ 
yses  in  the  absence  of  input  data. 

For  acoustic  freguencies  below  1000Hz  and  based  on 
results  from  the  FACT9H  Transmission  Loss  model,  the  errors 
in  the  EOTS  and  TEOTS  analyses  were  insignificant  for  (a) 
long  range  (convergence  zone)  predictions,  (b)  direct  range 
predictions  for  sonar  systems  with  a  figure  of  merit  of  90db 
or  greater,  and  (c)  direct  range  predictions  for  cross  layer 
and  below  layer  source/receiver  dispositions.  The  errors  in 
the  model  analyses  were  found  to  be  significant  for  sonar 
systems  with  a  figure  of  merit  of  80 db  or  less  when  both 
source  and  receiver  were  shallow. 

The  errors  in  modeled  mixed  layer  depth  were  found  to 
have  the  greatest  effect  on  acoustic  prediction  errors,  and 
the  effect  was  found  to  be  directly  proportional  to  acoustic 
freguency. 

A  comparison  between  the  FACT9H  and  PE  Transmission  Loss 
models  revealed  significant  differences  between  the  outputs 
of  the  two  models. 

The  PE  model  appeared  far  more  sensitive  to  errors  in 
the  temperature  profile  than  did  the  FACT9H  model.  The 
direct  range  errors  obtained  from  the  PE  model  outputs, 
based  on  modeled  and  measured  temperature  profiles,  showed 
no  dependence  on  the  figure  of  merit  and  less  dependence  on 
the  source/receiver  disposition  than  did  those  based  on  the 
FACT9B  model  outputs. 

The  disparity  between  FACT9H  and  ?E  model  sensitivities 
made  an  absolute  determination  of  the  operational 


significance  of  errors  in  the  EOTS  and  TEOTS  analyses 
inpossible.  However,  the  naxinun  effect  of  depth/teaperature 
profile  differences  on  direct  range  predictions  was 
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The  staff  at  FNOC  are  presently  restructuring  the  ocean 
thermal  analysis  aodel  algorithas  to  enable  the  aodels  to 
respond  aore  readily  to  input  data,  particularly  the  PLD 
segaent  cf  the  analysis.  The  physics  of  the  TEOTS  analysis 
are  also  being  incorporated  into  the  finer  grid,  regional 
analyses.  The  effectiveness  of  such  changes  should  be  tested 
in  a  siailar  areal  Banner  to  that  of  the  present  study,  and 
for  aany  aore  areas  of  the  world's  ocean,  to  establish  the 
operational  value  of  such  an  analysis  systea.  The  ainiaua 
data  input  reguired  to  give  a  representative  analysis  under 
varying  conditions  should  also  be  established,  in  order  that 
assets  aay  be  deployed  effectively  and  econoaically. 

For  regions  where  anoaalous  events  are  known  to  occur, 
the  feasability  of  utilizing  a  guasi-dynaaic  cliaatology, 
based  on  hindcasting,  should  be  studied,  i.e.,  climatology 
could  be  updated  when  a  persistent  anoaalous  event  is  known 
to  be  present.  This  would  give  a  aore  representative  anal¬ 
ysis,  at  depths  below  those  influenced  by  the  TOPS  elements 
of  the  TEOTS  analysis,  in  the  absence  of  input  data. 

Satellite  IB  iaagery  could  be  a  powerful  asset  in  data- 
sparse  areas  of  the  oceans,  providing  a  necessary  input  to 
ocean  theraal  analysis  aodels.  However,  the  results  of  the 
present  study  have  shown  that  the  relationship  between  SST 
and  sub-surface  features  is  not  always  an  obvious  one. 
Studies  should  be  undertaken  to  relate  surface  features  with 
sub-surface  structure,  either  through  direct  correlations  or 
aodeling.  Such  relationships  would  also  assist  the  planning 
of  saapling  strategies  aiaed  at  achieving  the  aost  econom¬ 
ical  and  efficient  use  of  assets  whilst  providing  sufficient 
input  to  the  ocean  theraal  analysis  aodel  to  achieve  the 
desired  accuracy. 
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The  effect  on  the  higher  frequency ;i.e. ,  greater  than 
iKhz,  active  sonar  performance  of  errors  in  the  model  anal¬ 
ysis  fields  should  be  examined.  The  loss  of  vertical  struc¬ 
ture  could  prove  important.  In  particular,  areas  where 
significant  vertical  features  are  present;  e.g.,  shallow 
sound  channels,  should  be  examined  and  the  feasibility  of  an 
additional  floating  level  in  the  analysis,  to  represent  the 
depth  of  such  features,  should  be  addressed. 

The  effect  of  mesoscale  variations,  on  scales  such  as 
those  observed  in  the  present  study,  on  the  other  parameters 
in  the  sonar  equation  should  be  examined;  e.g.,  ambient 
noise  changes  due  to  changes  in  the  marine  ecosystem  in  the 
region  of  the  weak  ocean  fronts  associated  with  upwelling 
regions. 

.  Operationally  important  areas  are,  in  general,  the  more 
densely  sampled  regicns  of  the  ocean.  However,  the  present 
fixed  regional  boundaries  utilized  in  the  FNOC  analysis 
models  are  only  broadly  representative  of  these  operational 
areas  and  generally  extend  beyond  the  limits  of  prime 
interest.  As  a  consequence,  the  shape  field  used  in  the 
analysis  is  influenced  by  possibly  less  well  sampled  areas 
outside  the  region  of  interest  with  possible  erroneous 
smoothing  effects  on  the  horizontal  and  vertical  thermal 
structure  that  the  analysis  would  otherwise  produce.  To 
remove  such  a  possibility,  the  utilization  of  a  flexible 
regional  boundary,  fitted  to  the  area  of  interest  ,  and  with 
a  grid  spacing  reprasentative  of  the  horizontal  variability 
of  the  region,  should  be  studied.  Such  an  approach  as  the 
one  outlined  above  would  also  facilitate  the  incorporation 
of  a  more  detailed,  possibly  dynamic,  climatology  into  the 
analysis. 
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